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Abstract 

A calculation of the expected signal due to Primakoff coherent conversion 
of solar axions into photons via Bragg scattering in several crystal detectors 
is presented. The results are confronted with the experimental sensitivities 
of present and future experiments concluding that the sensitivity of crystal 
detectors does not challenge the globular cluster limit on the axion-photon 
coupling g a ~ n - In particular, in the axion mass window m a > 0.03 eV explored 
with this technique (not accessible at present by other methods) <7 a77 might 
be constrained down to 10~ 9 GeV^ 1 (the recent helioseismological bound) 
provided that significant improvements in the parameters and performances 
of these detectors be achieved and large statistics accumulated. This bound 
should be considered as a minimal goal for the sensitivity of future crystal 
experiments. Consequently, finding a positive signal at this level of sensitivity 
would necessarily imply revisiting other more stringent astrophysical limits 
derived for the same range of m a values. 
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I. INTRODUCTION 



Introduced twenty years ago as the Nambu-Goldstone boson of the Peccei- 
Quinn symmetry to explain in an elegant way CP conservation in QCD 
the axion is remarkably also one of the best candidates to provide at least 
a fraction of the Dark Matter needed in Cosmology in order to explain both 
gravitational measurements and models of structure formation. 

Axion phenomenology is determined by its mass m a which in turn is fixed 
by the scale f a of the Peccei-Quinn symmetry breaking, m a ~ 0.62 eV (10 7 
GeV/f a ). No hint is provided by theory about where the f a scale should 
be. A combination of astrophysical and nuclear physics constraints, and the 
requirement that the axion relic abundance does not overdose the Universe, 
restricts the allowed range of viable axion masses into a relatively narrow 
window [§-f|]: 

10" 6 eV < ma < l(T 2 eV 
3 eV < m a < 20 eV. 

The physical process used in axion search experiments is the Primakoff effect. 
It makes use of the coupling between the axion field ip a and the electromag- 
netic tensor: 

£ = -^ga^ae^apF^F 013 = -g ail ^ a B ■ E 

and allows for the conversion of the axion into a photon. 

Like all the other axion couplings, g ail is proportional to m a f2],g||: 

9arn c 0.19^ (l - ^i±4) 10-GeV-i 
y 77 eV \N 3(1 + 2!)/ 

where E/N is the PQ symmetry anomaly and the second term in parenthesis 
is the chiral symmetry breaking correction. The anomaly E/N depends on the 
particular axion model, while the symmetry breaking correction is a function 
of the parameter z = mu/rrid ~ 0.56. Two popular models are the GUT- 
DFSZ axion § (E/N=8/3) and the KSVZ axion (E/N=0). However, it is 
possible to build viable axion models with different values of E /N ||] and the 
determination of the parameter z is subject to some theoretical uncertainties 
H . This implies that a very small or even vanishing g ail cannot be in principle 
excluded. 

Solid state detectors provide a simple mechanism for axion detection 



10 , 1 1 1 . Axions can pass in the proximity of the atomic nuclei of the crystal 
where the intense electric field can trigger their conversion into photons. In 
the process the energy of the outgoing photon is equal to that of the incoming 
axion. 

Axions can be efficiently produced in the interior of the Sun by Primakoff 
conversion of the blackbody photons in the fluctuating electric field of the 
plasma. The resulting flux has an outgoing average axion energy E a of about 



4 keV (corresponding to the temperature in the core of the Sun, T ~ 
that can produce detectable x-rays in a crystal detector. Depending on the 
direction of the incoming axion flux with respect to the planes of the crystal 
lattice, a coherent effect can be produced when the Bragg condition is fulfilled, 
leading so to a strong enhancement of the signal. A correlation of the expected 
count -rate with the position of the Sun in the sky is a distinctive signature of 
the axion which can be used, at the least, to improve the signal/background 
ratio. 

The process described above is independent on m a and so are the achiev- 
able bounds for the axion-photon coupling (7,177 • This fact is particularly ap- 
pealing, since other experimental techniques are limited to a more restricted 
mass range: "haloscopes" [12], that use electromagnetic cavities to look for 
the resonant conversion into microwaves of non relativistic cosmological dark 
halo axions, do not extend their search beyond m a ~ 50 fieV, while the dipole 
magnets used in "helioscope" [13| experiments are not sensitive at present to 
solar axions heavier than m a ~ 0.03 eV. 

A pilot experiment carried out by the SOL AX Collaboration 14] has al- 
ready searched for axion Primakoff conversion in a germanium crystal of 1 kg 
obtaining the limit g ail < 2.7 x 10~ 9 GeV -1 . This is the (mass independent 
but solar model dependent) most stringent laboratory bound for the axion- 
photon coupling obtained so far, although less restrictive than the globular 
cluster bound [15] g ail < 0.6 x 10~ 10 GeV -1 . Notice however that the exper- 
imental accuracy of solar observations is orders of magnitude better than for 
any other star. 

As a matter of fact, the solar model itself already requires [16| g a ~ n < 10~ 9 
GeV -1 , whereas the above mentioned Ge experiment has not yet reached such 
sensitivity. The 10 _9 GeV~ 1 limit sets, then, a minimal goal for the sensitiv- 
ity of future experiments, prompting the need for a systematic discussion of 
present efforts and future prospects for axion searches with crystals. In the 
following we give the result of such an analysis, focusing on Germanium, Te02 
and Nal detectors. 



II. PRIMAKOFF CONVERSION IN CRYSTALS 

We will make use of the calculation of the flux of solar axions of Ref. 
| |17| |, with no direct coupling between the axion and the leptons (hadronic 
axion), with the modifications introduced in Ref. [11] to include helium and 



metal diffusion in the solar model. A useful parametrization of the flux is the 
following: 

d$_ f-Qp (E a /E f 

dE a VA E e E «/ E °-l {) 

where A=(<7 a77 x 10 8 /GeV -1 ) 4 is a dimensionless coupling introduced for later 
convenience, &q=5.95 xlO 14 cm -2 sec -1 and I?o=l-103 keV. 
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The expected count-rate in a solid-state ionization detector, integrated 
in the electron-equivalent energy window Ei<E ee <E2, is calculated starting 
from the expression: 



pE 2 roo 

R(E h E 2 ) = dE ee dE 1 -—{E 
Je 1 Jo dE^ 



rE 2 rco i (B ee -fi 7 ) 2 



7 AV27T 



where EL is the energy of the outcoming photon, A is the resolution of the 
detector, FWHM ~ 2.35 A, while: 

a is the cross section of the process. The recoil of the nucleus can be neglected, 
so the energy of the outgoing photon is equal to that of the incoming axion 
and the differential cross section for the Primakoff conversion may be written 
as: 



dELdft dtt 

where [18 



29 and is the scattering angle, while 

F a (q) = k 2 J d^x^xY** 

is the Fourier transform of the electrostatic field 4>. q is the transferred mo- 
mentum, q = \q\ = 2ksin9 and k=\k\^E a is the axion momentum. 

The energy distribution of Eq.(||) implies that q corresponds to a wave- 
length of a few A, which is of the order of the distances between atoms in 
a typical crystal lattice. As a consequence of this a Bragg-reflection pattern 
arises in the calculation of the Fourier transform of Eq.(||) due to the periodic 
properties of the electric field in the crystal. Using translational invariance 
the periodic field may be written as [O] (we consider here the general case of 
a multi-target species): 

ij ij 1 *l G 

where Xi are the positions in space of the elements of the lattice, Zj is the 
atomic number of the j— th type target nucleus, rj is the screening length 
of the corresponding atomic electric field parametrized with a Yukawa-type 
potential, and the sum in the last term is intended over the vectors G of the 
reciprocal lattice, defined by the property expzGj • = 1. A multi-target 
crystal is described by a Bragg lattice with a basis whose sites are occupied 
by atoms of different types. Indicating by a\ the i'th basis vector occupied 



by the j'th target-nucleus type, the uq coefficients appearing in Eq.(lO) turn 
out to be: 



where 

i 

are the structure functions of the crystal and 



Zjek 2 
r~ 2 + q 2 



is the form factor of Eq.(||) for a single j-type target. v a indicates the volume 
of the elementary cell of the lattice. 

Integration of Eqs.(j5|,P) over E ee and E a leads to: 



R(E 1 ,E 2 ) = / 2c 



d 3 q d$ 



( Ea )^\F a (q)\ 2 sin\29)W 



16vr 2 ' 



CJ77 



t<Z dE a 



|GP16 

where V is the volume of the detector, 



Y J Fl j {G)S 3 {G)\ 2 W{26)W 



W(E u E 2 ,E a ,A) 



1 



erf 



Eg - Ei 
V2A 



erf 



Eg — Ei 
V2A 



and in Eq.(14) |.F a (<7)| 2 has been expanded by making use of Eqs.(p|,10|,ll,l^). 
In the final result the integral over the transferred momentum has been re- 
placed by a sum over the vectors of the reciprocal lattice, i.e. over the peaks 
that are produced when the Primakoff conversion verifies the Bragg condition 
q = G and the crystal interacts in a coherent way. 

The Bragg condition implies that in Eq.(|l5|) E a = ^ c ^^ where the uni- 
tary vector u points toward the Sun. This term induces a time dependence 
in the expected signal. 



III. TIME CORRELATION AND BACKGROUND REJECTION 

In the signal the dependence on A can be factorized: R = XR. An example 
of the function R for several materials is shown in Fig. 1 as a function of time 
during one day and assuming the coordinates of the LNGS laboratory. The 
crystallographic inputs in the calculation are given in the Appendix. In all 
the plots the main axis of the lattice have been taken parallel to the South, 
West and upward directions. This explains the symmetry shown by all the 
plots. 



The signal is peaked in energy around the maximum of the flux of Eq. (Q) 
and presents a strong sub-diary dependence on time, due to the motion of 
the Sun in the sky. The time duration of the peaks decreases with growing 
energies, from tens of minutes in the lowest part of the axion energy window, 
down to a minimum of about one minute in the higher one, and is related to 
the energy resolution of each detector. 

In order to extract the signal from the background for each energy interval 
Ef. < E < Ef. + AE we introduce, following Ref. [O], the quantity: 



n 



X = E <R>]-n i = J2W i -n i 

i=l i=l 

where the indicate the number of measured events in the time bin ti,ti + At 
and the sum is over the total period T of data taking. Here and in the following 
the brackets indicate time average. 

By definition the quantity x is expected to be compatible with zero in ab- 
sence of a signal, while it weights positively the events recorded in coincidence 
with the expected peaks. 

The time distribution of n« is supposed to be given by a Poissonian with 
average: 

< m >= [XR(ti) + b] At. 

Assuming that the background b dominates over the signal the expected av- 
erage and variance of x are given by: 

<X> = X-A 
a 2 ( X )^b/A 

with A = J2i WfAt. Each energy bin E k ,E k + AE with background b k 
provides an independent estimate A^ = Xk/Ak so that one can get the most 
probable combined value of A: 



a(X)= \^2A k /b k j . 

The sensitivity of an axion experimental search can be expressed as the 
upper bound of g a ^ which such experiment would provide from the non- 
appearance of the axion signal, for a given crystal, background and exposition. 
If A is compatible to zero, then at the 95% C.L. A < 2x 1.64xcr(A). 

It is easy to see that the ensuing limit on the axion-photon coupling g£™ 
scales with the background b (assumed to be flat) and exposure MT in the 
following way: 
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where M is the detector mass and T the total time. The factor K~ 6.5 
[Efc < (W(E k )/M) 2 > keV/AE k ]~ l/8 is a function of the parameters of the 
crystal, as well as of the experimental threshold and resolution. The ap- 
plication of the statistical analysis described above results in a background 
rejection of about two orders of magnitude. 



IV. DISCUSSION AND CONCLUSIONS 

In order to perform a systematic analysis of the axion-detection capability 
of crystal detectors, we have applied the technique described in the previous 
section to several materials. The result is summarized in Table |, where the 
limit given by the experiment of Ref. [14] is compared to those attainable 



with running |2(|, being installed [21 ,22] and planned [p5|,24| crystal detector 
experiments. 

In Table | a Pb detector is also included, to give an indication of the 
best improvement that one would expect by selecting heavy materials to take 
advantage of the proportionality to Z 2 of the cross section of Eq.(|8|). These 
results can be easily extended to other elements with lower atomic numbers 
than those considered in Table 1 (for instance in AI2O3 or LiF crystals, also 
in operation or planned dark matter experiments) although they are expected 
to yield less stringent limits. 

In Fig. 2 the result of our analysis is compared to the present astrophysical 
and experimental bounds in the plane m a -g a ^. The horizontal thick line 
represents the constraint g a ^ <3xl(r 10 GeV" 1 taken from Table | The 
mass intervals m a < 10~ 6 eV, 1CP 2 eV < m a < 3 eV, m a > 20 eV are 
excluded respectively by cosmological overclosure Q, SN1987A || and oxygen 
excitation in water Cherenkov detectors @, the edges of all the excluded 
regions being subject to many astrophysical and theoretical uncertainties. 
The horizontal lines represent the limits from SOLAX (T^[, helioseismology 
[16], and globular clusters |15[ ] . The theoretical predictions for g a77 in the 



case of the DFSZ and the KSFZ axion models are also given. 

As shown in the expression of the g a77 bound of Eq.(p2[), the improvement 
in background and accumulation of statistics is washed out by the 1/8 power 
dependence of <7 a77 on such parameters. It is evident, then, that crystals 
have no realistic chances to challenge the globular cluster limit. A discovery 
of the axion by this technique would presumably imply either a systematic 
effect in the stellar-count observations in globular clusters or a substantial 
change in the theoretical models that describe the late-stage evolution of 
low-metallicity stars. 

On the other hand, the sensitivity required for crystal-detectors in order 
to explore a range of g a -yy, compatible with the solar limit of Ref. [16], ap- 



pears to be within reach, although only marginally, and provided that large 
improvements in background as well as substancial increase of statistics be 
guaranteed. 
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In light of the previous discussion, prospects to detect the axion with 
crystal detectors appear to be dim. The only realistic possibility to exploit 
this technique seems to rely on the fact that collecting a statistics of the order 
of a few tons x year could be achievable by adding properly the results of the 
various dark matter search experiments already planned for the near future. 
As shown in Fig. 2 the ensuing limit for crystals would cross the theoretical 
models for an axion mass in the range of a few eV, in a region compatible 
with the bound given by SN1987A. 



APPENDIX 

In this section we summarize the crystallographic information necessary 
to calculate the axion-photon coherent conversion discussed in section || for 
the materials shown in Table pL We recall that a Bravais lattice generated 
by the primitive vectors a?, has a reciprocal lattice generated by the basis 
Gi = 2ir/v a €ijkXj x Xk, where v a = x\- (a?2 x £3) is the volume of the primitive 
cell, while e^fe is the totally anti-commutating tensor. 

Germanium, sodium iodide and lead share the same type of underlying 
Bravais lattice, the faced-centered cubic crystal (fee) structure generated by 
the primitive vectors 

xi = |(0,l,l) 

X 2 = |(1,0,1) 
X3 = |(l,l,0). 

The corresponding primitive cell sizes a are given in Table |ll[ along with the 
basis vectors for Ge and Nal. 

The Bravais lattice of Te02 is tetragonal. In the plots of Fig. 2 the longer 
side of the cell has been chosen in the vertical direction. In this case the 
generators are simply: 

xi = a(l,0,0) 
x 2 = o(0, 1,0) 
£3 = 6(0,0,1). 

with a=4.796 A and b=7.626 A. Each primitive cell contains a basis of 4 
Tellurium atoms and 8 Oxygen atoms. Their positions in the primitive cell 



are given in Table 111 



The parameter r*j entering the Yukawa potential of Eq.(10) is of the order 
of the atomic dimensions. In our calculation we have assumed rj=l A for all 
the crystals. Varying rj in the range 0.5 A < r j < 1 A affects the g a ^ limit 
by less than ~ 30%. 
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TABLES 



TABLE I. Axion search sensitivities for running (DAMA f20|), being installed (CUORICINO 
1 21|, CANFRANC ||) and planned (CUORE flEJ, GENIUS @]) experiments are compared to 
the result of SOLAX fll4| . To estimate the best improvement that would be possible by selecting 
high-Z materials, a calculation is also shown for a Pb lattice, assuming the same conditions of 
background, threshold and resolution of the GENIUS experiment. The coefficient K has been 



introduced in Eq.(22). 





K 


M 
(kg) 


b 

(cpd/kg/keV) 


Eth 
(keV) 


FWHM 

(keV) 


g l j™(2 years) 
(GeV- 1 ) 


Ge 


14 


1 


2.5 


1 


3 


4 


1 


2.7xl0~ 9 


Ge [ 


24 


1 


2.5 


1000 


lxlO" 4 


4 


1 


3xl0" 10 


Te0 2 


m 


3 


42 


0.1 


5 


2 


1.3xl0" 9 


Te0 2 


1 


I 


2.8 


765 


1 x 10~ 2 


3 


2 


6.3xl0" 10 


Nal 


11 


2.7 


87 


1 


2 


2 


1.4xl0" 9 


Nal 


© 


2.8 


107 


2 


4 


2 


1.6xl0" 9 


Pb 


2.1 


1000 


lxlO" 4 


4 


1 


2.5xl0" 10 



TABLE II. Crystallographic information for Ge, Nal, and Pb: a is the primitive cell size and 
the a\ are the basis vectors entering in the evaluation of the structure functions of Eq.(12). 



Species 


a (A) 


Basis a\ 


Ge 


5.66 


aj = (0,0,0) 

4=f(i,i,i) 


Nal 


6.47 


a}=(0,0,0) (I) 
3?=§(1,1,1) (Na) 


Pb 


4.95 


No basis 



TABLE III. Basis vectors for Te0 2 entering in the evaluation of the structure function of 



Eq.(12). The longer side of the tetragonal primitive cell has been chosen in the vertical direction. 
All numbers are in A. 



Oxygen 


Tellurium 


a\=( 0.998 , 0.916 , 1.480) a\={ 3.797 , 3.879 , 5.293) 
4=( 1-481 , 3.396 , 3.386) a\={ 3.314 , 1.399 , 7.199) 
a\={ 1.399 , 3.314 , 0.426) a\={ 3.396 , 1.481 , 4.239) 
4=( 0.916 , 0.998 , 6.145) a\=\ 3.879 , 3.797 , 2.332) 


a\={ 0.093 , 0.093 , 0.000) 
al=( 4.702 , 4.702 , 3.813) 
a§=( 2.304 , 2.491 , 1.906) 
d\={ 2.491 , 2.304 , 5.719) 
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FIGURE CAPTIONS 



Figure 1 - Expected axion signals for Primakoff conversion in various 
crystals as a function of time for A = 1. In the calculation the representative 
day of 1 april 1998 and the coordinates of the LNGS laboratory have been 
assumed. From top-left to bottom-right: a) Ge, 2 keV< E ee <2.5 keV; b) Ge, 
4 keV< E ee <4.5 keV; c) Te0 2 , 5 keV< E ee <7 keV; d) Te0 2 , 7 keV< E ee <9 
keV; e) Nal, 2 keV< E ee <4 keV; f) Nal, 4 keV< E ee <6 keV. 

Figure 2 - The solar axion limit attainable with crystal detectors (hori- 
zontal thick line) is compared to the present astrophysical and experimental 
bounds and to the DFSZ and KSVZ axion theoretical predictions. The shaded 
mass intervals m a < 10 -6 eV, 10~ 2 eV < m a < 3 eV, m a > 20 eV are ex- 
cluded respectively by cosmological overclosure ||, SN1987A Q and oxygen 
excitation in water Cherenkov detectors ||]. The horizontal lines represent the 
upper limits for g a -y^ from helioseismology Q (dashes), SOLAX Q (dots) 
and globular clusters [|D| (dot-dashes). The regions excluded by haloscope 
searches [12| and the upper limit on g a ~ n from the Tokio Helioscope experi- 
ment [13 1 are also shown. 
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